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-  INTRODUCL^ZQK" 

The  existing  IBM  7C90  computer  program  to  calculate  static 
stress  distributions  of  propeller  blades  (ref. 2)  has  been 
revised  to  execute  the  full  potential  of  the  program,  i.e. 
to  generate  108  degrees  of  freedom.  As  noted  in  the  con¬ 
clusions  of  refo  2,  the  round-off  error  accumulated  in  the 
process  of  solving  a  system  of  more  than  90  equations  had 
precluded  the  full  exploitation  of  the  existing  program. 

The  aforementioned  revision  consisted  of  developing  a 
technique  which  voiild  result  in  a  solution  for  a  system 
as  high  as  108  equations. 

Two  techniques  were  investigated.  One  was  to  utilize 
a  ’’brute  force"  approach  and  invert  the  coefficient  matrix 
in  double  precision,  using  a  routine  which  ignored  the  fact 
that  the  matrix  was  positive-definite.  The  other  approach 
v/as  to  pai'tition  the  matrix  into  9  sub-matrices  (6  independ¬ 
ent  due  to  siTTimet-ry)  and  generate  the  solution  vector; 
the  inverse  is  not  calculated.  Either  technique  would 
then  be  followed  by  appropi  .ate  iteration  schemes  which 
■v.'Ould  refine  the  solution  vector  to  the  required  precision. 

ThiO  direct  inversion  method  enjoyed  the  advantage  of 
alv7ays  being  able  to  generate  a  solution  vector;  past 
experience,  however,  indicated  that  direct  inversion  of 
such  a  large  system  invariably  produced  large  round-off 
errors  and  a  corresponding  relatively  inaccurate  solution 
vector.  It  was  necessai’y  to  determine  whether  the  sub¬ 
sequent:  iteration  schemes  would  filter  out  the  "noise"  in 
this  solution  within  a  reasonable  amount  of  machine  time. 
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':che  Gauss-Seidel  iteration  is  unconditionally  convergent  for 
a  positive-definite  matrix;  however,  convergence  can  be 
extremely  slow  -  perhaps  prohibitively  slow  -  if  the  trial 
solution  vector  is  a  very  poor  approximation  of  the  true 
sciuticn. 

The  partitioning  method  operates  on  a  family  of  sub¬ 
matrices  of  size  36  x  36  at  most,  rather  than  a  single 
108  X  108  matrix.  Three  inversions  are  executed  on  these 
smaller-sized  matrices;the  resulting  round-off  error  using 
do'able-'precision  arithmetic  is  virtually  nill.  However, 
as  shown  in  detail  later,  the  subsequent  operations  include 
subtractions  of  matrices  of  like  terms;  it  is  here  that  the 
introduction  of  noise  is  likely  to  occur.  Thus,  while  a 
solution  vector  is  always  generated,  its  accuracy  is  sub¬ 
ject  to  question. 

This  report  should  be  viewed  as  a  modification  to  ref. 2. 
All  symbols  are.  defined  in  ref.  2  and  descriptive  figubes 
arc  contained  therein. 
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II  ANALYSIS 

A-  Direct  Inversion  Method 

Given  the  system  Ax  =  f,  the  inverse  of  the  co¬ 
efficient  matrix,  A“^  is  calculated  leading  to  the 
solution  vector,  x  =  A'^  f.  The  inverse  actually  gen¬ 
erated  is  in  error  due  to  round-off  error;  call  it  A”^. 
Thus,  the  resulting  solution  x  =  A"^  f  is  also  in  error. 

If  wc  write  x  =  x  +  Ax  and  define  f  =  Ax  and  f  =  ?  +  A  f , 
then  we  may  write  the  system  as 


A  (  X  +  A  x)  =  (f  +  A  f) 

which  leads  to  A  A  x  =  A  f,  where  A  f  =  f  -  f  may  be  calcula- 

^  -1 

ted.  If  the  inverse,  A  is  "sufficiently"  accurate  (ref. 3), 
then  we  may  set  up  an  iterative  scheme  as  follows: 


A  A  f 


^(n+1)  ^  ^(n)  ^  ^  ^(n)_  ^(n+1)  ^  ^  f  .  f(n+l) 


Thus,  if  the  iteration  is  convergent.  Ax  approaches  zero 
and  the  desired  precision  of  the  solution  vector  is  attained. 
Subsequently,  a  Sidelian  iteration  is  performed  to  further 
refine  the  solution.  During  the  course  of  each  sweep  thru  the 
systen),  the  elements  of  the  solution  vector  are  calculated 
in  terms  of  the  "most  recent"  values: 


x,  = 


aij  Xj 


1  —  N 


y=i 
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In  the  above  discussion^ 

A  is  the  matrix  of  coefficients  (strain  energy 

matrix,  ref.  1) 

X  is  the  (unknovvn)  solution  vector 

f  is  the  forcing  vector  which  is  calculated 

previously  by  the  program,  utilizing  the 
pressure  loading  input  data. 

N  is  the  size  of  the  system  ^108 

Ax,  A  f  are  error  vectors 


B.  Partitioning  Method 


The  system  Ax  =  f  may  be  partitioned  into  sub-matrices, 
^ij,  and  vectors,  Xj^  and  f^  as  follows: 


These  partitioned  elements  may  be  operated  upon  by  a 
Causs-Jordan  reduction  scheme  to  yield  an  upper  triangular 
system; 


r 

"7 

i 

f  1 

i 

r  1 

^11 

^12 

®13  i 

j 

^1 

i^i 
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where 


^22 

®22 

'  ^21 

-u 

^12 

^23 

"  ^23 

“  ^21 

-1 

^11 

^13  ? 

T 

^32  -  ^23 

^33 

=  ^33 

“  ^31 

-1 

^11 

ai3 

=33 

”  ^33 

”  ^32 

^22“^ 

^23 

^2 

=  f2 

^21 

-1 

^11 

fl 

^3 

=  ^3 

"  ^31 

-1 

^11 

fl  - 

-1 

h32  b22  ^2 

Then  utilizing  "back  substitution" , 

“•1 

^3  ■  =33  ^3 

-1  r  '1 

=  ^22  '^2  “  ^23 

i*  ^ 

r  ■ 

xi  =  a^i  -ai2  ^2  -  ai3X3  .+ 


Note  that  since  the  coefficient  matrix.  A,  is  symmetric, 

T 

is  symmetric  and  c\jj^  =  aj_j  i  There  are  3  inversions 
performed  -  submatrices,  ^22*  *^33  “  inverse 


-1 


is  not  generated. 


After  the  solution  vector  is  calculated  the  program 
utilizes  rhe  Sidelian  iteration  scheme  to  refine  the 
results,  as  described  previously. 
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III  RESULTS  OF  TEST  RUNS 

Two  program  decks  were  formed  and  were  tested  in 
parallel.  Using  the  DD828  propeller  blade  configuration 
as  the  test  vehicle,  these  programs  were  run  on  the 
DTMB  IBM  7090  computer.  To  study  the  influence  of  error 
growth  as  the  size  of  the  system  was  increased,  the 
programs  utilized  30,  45,  63  then  84  degrees  of  freedom, 
in  that  sequence.  The  results  of  the  two  versions  were 
compared,  and  these  in  turn  were  compared  with  the  results 
obtained  with  the  previous  version  of  the  program. 

All  results  were  in  excellent  agreement  with  the  4.0ne 
exception  of  the  partitioning  version  at  84  degrees  of 
freedom.  The  two  new  programs  were  then  run  at  75  degrees 
of  freedom,  and  it  was  found  that  several  elements  of  the 
solution  vectors  of  the  two  new  versions,  differed  signif¬ 
icantly.  Further  examination  indicated  conclusively  that 
the  partitioning  version  was  not  reliable  beyond  63  degrees 
of  freedom. 

The  direct  inversion  program  was  then  run  at  full 
capacity  -  108  degrees  of  freedom  -  and  the  results  com¬ 
pared  favorably  with  those  of  84  degrees  of  freedom.  To 
est  the  accuracy  of  solution,  the  number  of  iterations 
utilized  by  the  Sidelian  technique  was  doubled.  This 
final  run  indicated  that  doubling  the  number  of  iterations 
did  not  alter  the  solution. 

These  results  have  been  plotted  and  appear  on  the 
following  pages. 


O'Z 
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IV  CONCLUSIONS 

The  computer  program  now  has  the  capability  of 
utilizing  108  degrees  of  freedom  in  the  process  of 
generating  stress  distributions  for  propeller  blades. 

The  results  of  such  a  test  run  have  been  plotted  and  may 
be  compared  with  similar  plots  of  results  which  were 
generated  using  84  degrees  of  freedom  (see  ref.  2) 

These  comparisons  indicate  that  no  radical  change 
in  stress  contours  occurs  as  a  result  of  i  icreasing  the 
number  of  degrees  of  freedom  from  84  to  108?  indeed,  there 
is  a  marked  similarity  between  the  two  sets  of  plots.  This 
similarity  indicates  that  the  solution  has  'Settled  down” 
and  that  108  degrees  of  freedom  are  adequate  for  the 
purpose  of  depicting  the  variation  of  displacement 
components  over  the  planform  of  propeller  blades.  Thus, 
increased  confidence  is  warranted  in  the  capability  of 
the  subject  computer  program  to  generate  accurate  stress 
distributions  for  static,  distributed  loadings  on  propeller 
blades. 


V  INPUT  FORMAT 
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Card 

1 

2 

3 

4 ,  7 /  10 ,  • « • 

5 ,  8 ,  11 , . . . 

6,  9,  12,... 

T  -  1 

T 

T+1,  T+3,  T+5,... 


T+2,  T+4,  T+6,... 


C 


Input  Information 
Identification 
Blank 

Poisson's  Ratio,  Blade  Root  Geometry 

Section  Properties 

Offsets  to  Back  Face  Along  Section 

Offsets  to  Front  Face  Along  Section 

Properties  of  Blade  Tip 

Code  Card 

Pressure- jumps  at  points  on  sections 
corresponding  respectively  to  those 
sections  of  input  cards  3,  6,  9,  ... 
These  points  are  from  leading  edge 
to  25%  of  chord. 

Pressure-jumps  at  points  on  sections 
corresponding  respectively  to  those 
sections  of  input  cards  3,  6,  9,  ... 
These  points  are  from  30%  of  chord 
to  trailing  edge. 

Control  Card 


Next  set  of  inputs,  if  any. 

Card  No.  1; 

Any  identification  state-  ent  of  71  characters  or  less  will 
be  accepted.  The  statement  must  be  punched  in  columns  2-72  of 

a  card.  The  digit,  "1"  should  be  punched  in  column  1  for 
convenience  of  output. 

Card  No.  2: 


Blank 
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Card  No.  3 


Data 

Columns 

Format 

V 

1-5 

X.XXX 

^R 

7-12 

XX. XXX 

14  -  19 

XX. XXX 

21  -  26 

XX. XXX 

28  -  33 

XX. XXX 

N 

37  -  38 

XX 

The  first  word,  v,  is  Poisson's  Ratio.  The  next  4 

words,  are  indicated  in  Pig.  2  (ref.  2). 

number  of  degrees  of  freedom  for  each  displacement  com¬ 
ponent  is  N.  Thus  for  63  degrees  of  freedom,  N  is  21. 

If  a  full-scale  problem  is  being  run  (108  degrees  of 
freedom) ,  N  may  be  left  blank. 

Card  No.  4 

Data  Columns  Format 


%  Radius 

1  - 

2 

XX 

Radius 

4  - 

9 

XX. XXX  or 

XXX. XX 

Span 

11  - 

16 

XX. XXX  or 

XXX. XX 

Pitch 

18  - 

24 

+XX.XXX  or 

+XXX.XX 

Rake 

26  - 

32 

+X.XXXX 

Skew 

34  - 

40 

lX.XXXX 

The 
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Card  !^o.  5 
Data 
A 
E 
F 
G 
H 
J 
K 
L 
M 
N 
P 

Card  No,  6 

Identical  in 

Card  T-1 
Data 


Radius 

Thickness 

Pitch 

Rake 

Skew 

Card  T 
Data 


Columns 
1-6 
7-12 
13  -  18 
19  -  24 

25  -  30 
31  -  36 
37  -  42 
43  -  48 
49  -  54 
55  -  60 
61  -  66 

format  to  Card  No, 

Columns 
1-2 
4-9 
11  -  16 
18  -  24 

26  -  32 
34  -  40 

Column 
1 


Format 

iX.XXX  or  +XX.XX 


Format 

Must  be  blank 
XX.X3DC  or  XXX. XX 
XX. XXX 
XXX. XXX 
+X.XXXX 
iX.XXXX 

Format 


Code 
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Card  T  +  1 

Data 

Columns 

Po 

1-7 

^1.25 

9-15 

P2.5 

17  -  23 

P5.O 

25  -  31 

^7.5 

33  -  39 

^10 

41  -  47 

Pl5 

49  -  55 

^20 

57  -  63 

^25 

65  -  71 

Card  T  +  2 

Data 

Columns 

P 

30 

1-7 

^40 

9-15 

^50 

17  -  23 

^60 

25  -  31 

^70 

33  -  39 

CO 

0 

41  -  47 

P90 

49  -  55 

^95 

57  -  63 

PlOO 

65  -  71 

Format 


XXX. XXX 


t 


J 


Format 

XXX.XXX 


Card  C 
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Data  Column  Format 

Code  1  X 

The  program  can  accept  a  maximum  of  12  sections  cor¬ 
responding  to  a  total  of  66  input  cards  per  set. 

Cards  5  and  6  specify  the  thickness  offsets  along  a 
section  at  equal  chord  increments,  from  leading-edge  to 

tr ailing-edge.  The  offset  of  A  is  the  offset  near  the  lead¬ 
ing  edge,  the  offset  of  E  is  at  one-tenth  the  chord  distance 
from  the  leading  edge,  etc.  The  offset  of  P  refers  to  trailing 
edge. 

Card  T  specifies  the  type  of  peaked  loading  if  any.  The 
code  is  *'l"  for  a  "peaked"  load  and  a  fully  wetted  blade,  and 
blank  for  any  load  which  is  not  "peaked". 

Cards  T+1  and  T+2  specify  the  complete  pressure-jump 
loading  acting  on  that  section  described  on  cards  4,  5  and  6. 

The  subscripts  of  the  p's  are  the  percentages  of  chord,  count¬ 
ing  from  leading  edge  to  trailing  edge.  For  example,  Pq  is  the 
pressure-jump  loading  acting  at  the  leading  edge  of  the  section 
under  consideration;  P20  acts  at  20%  of  chord,  or  one-fifth  the 
distance  from  the  leading  edge.  Of  course,  is  the  pressure 

jump  loading  at  the  trailing  edge  (invariably,  zero) .  Cards 
T+3  and  T+4  specify  the  pressure- jump  loading  on  that  section 
described  on  cards  7,  8,  and  9;  cards  T+5  and  T+6  pertain  to- 
cards  10,  11,  and  12,  etc. 
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The  control  card  dictates  the  action  of  the  program 
after  a  "run"  has  been  completed.  If  no  more  inputs  follow, 
the  control  card  may  be  blank.  If  it  is  desired  to  conduct 
an  additional  analysis  on  the  same  blade  using  a  different 
loading,  the  control  Ccird  must  have  a  "2"  punched  in  Column  1, 
and  be  followed  immediately  by  cards  specifying  the  loading, 
as  described  above  (cards  T  through  C) .  Should  it  be 
desired  to  conduct  an  analysis  on  another  blade  (no  matter 
what  the  loading  may  be),  a  "1"  must  be  punched  in  Column  1, 
followed  by  cards  describing  the  new  geometry  and  correspond¬ 
ing  pres sure- jump  loading  (cards  1,  2,  3...,C). 

The  form.ats  indicated  are  representative.  The  decimal 
points  may  be  shifted  in  either  direction,  so  long  as  the 
input  "word"  is  confined  within  the  columns  specified. 

The  program  assumes  that  the  tip  of  the  blade  is  a 
point  (span  equal  to  zero)  as  is  the  case  of  a  propeller 
blade.  If  the  section  at  the  blade  tip  has  some  span, 
such  as  the  flat  plates  it  will  be  necessary  to  input 
two  sections  extremely  close  to  the  tip.  In  general,  both 
sections  should  be  within  the  outer  1%  of  the  total  blade 


span. 
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